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Abstract 



The aim of this work is to show how symbolic computation can be used to perform 
multivariate Lagrange, Hermite and Birkhoff interpolation and help us to build more 
realistic interpolating functions. After a theoretical introduction in which we analyze the 
complexity of the method we shall focus our attention on applications. 



Introduction 

Multivariate interpolation consists in given finitely many points pi, . . . ,pt in a n-dimensional 
affine space, over a field K, such that pi / pj if i / j, and t scalars vi, . . . ,vt, determining 
a polynomial F £ K[Xi, . . . ,X„,] such that F{pi) = vi for any i = 1, . . . ,t. This is known 
as the Lagrange interpolation problem. It is well known that there are several methods and 
algorithms to compute such a polynomial F. 

In this paper we will add some extra conditions and show some algorithms to find an interpo- 
lating polynomial. We will deal with possible values of the interpolating polynomials at the 
given points, as in the Lagrange problem, and also with values of some higher order derivatives 
at these points. We consider two different cases. The first deals with considering unitary vec- 
tors U in iiT" and values for all derivatives DjjF{pi) for j = 1, . . . , s{i), the so called Hermite 
interpolation problem consists in determining a polynomial F G K[Xi, . . . ,Xn] satisfying all 
these conditions, see OS] and [6]. The second deals with the same situation, but we assume 
some gaps may exist, i.e., there may exists some point pi, some order j < Si or some vector U 
such that we have no assigned value for DjjF{pi); the so called Birkhoff interpolation problem 
consists in determining a polynomial F £ K[Xi, . . . ,Xn] satisfying these conditions. 

In two of the three cases when we homogenize" the relationships, i.e., when we consider 
^jjiPi) = ^iid F{pi) = for any U, j and pi, the set of all polynomials satisfying all 
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these relationships constitutes an ideal of the polynomial ring K[Xi, . . . This allows 

us to use some computational method in K[Xi, . . . at least when K is a computable 

field, in order to determine a general solution to the interpolation problem. This method is 
based on the application of Groebner basis. On this matter we study some results about the 
complexity of the algorithm and the best monomial order we need to use. The third case, 
the Birkhoff interpolation problem, does not produce ideals in the polynomial ring, hence 
the application of Groebner basis is not allowed. But we may reorganize the information to 
work as in the previous two cases to get a general solution of the interpolation problem, and 
after that, to particularize to get an specific one. Nevertheless in this case every gap will 
introduce a parameter to the specific solution, hence we will obtain the Birkhoff interpolation 
polynomial as a solution to an indeterminate compatible system of linear equations. In the 
practice, a high number of parameters (or equivalently gaps) slows down the computation of 
the interpolation polynomial. 

One of the theoretical problems we solve in this paper is to determine the existence and 
uniqueness of the Hermite interpolation polynomial, which in the Lagrange case is well known. 
We obtain this result as a consequence of the study of solutions to the Birkhoff interpolation 
problem. 

Let us describe the content of each section in this paper. 

In Section one we recall some facts about multivariate Lagrange interpolation, and in Section 
two we introduce the multivariate Hermite interpolation by using Groebner basis theory. In 
particular we prove that if we consider elements pi, ... ,pt in A'^{K) such that pi / pj if i ^ j, 
unitary vectors Uij- G K'^, i = 1 . . . , t, jj = 1 . . . , s{i), and sets Hi C W^^^ such that if h £ H 
and /i — e G then h — e G Hi, and define for any h = {hi, . . . ,hg(^i)) G Hi and any 

polynomial F G K[Xi, ...,Xn] the derivative D^^^F = Z)[|^\^ • • • D^I^^^'^^F, then for any set 

{vi^h G -ftT I /i G Hi} there is a polynomial F G K[Xi, . . . ,Xn] satisfying D^^^F{pi) = Vi^h 
for any pi and any h G H^. This polynomial F is unique modulo the cofinite ideal {F G 
K[Xi, ...,Xn]\ D^^)F{pi) = 0, for all pi and he Hi}. 

In Section three we consider the more intricate case in which some gaps appear, i.e., when 
the sets Hi have some gaps. In general the set {F G K[Xi,...,Xn] \ D^'''^F{pi) = 
0, for all Pi and h £ Hi} is not an ideal, but we may extend the set of conditions adding 
those corresponding to the gaps. In this case we have new sets. Hi, in such a way that 
the set {F G K[Xi, . . . ,Xn] \ D'^'^'>F{pi) = 0, for all and h G Hi} is a cofinite ideal. 
Let G be a Groebner basis, then a K-Yectoi space basis of the quotient ring is parameter- 
ized by N" \ Exp(G), and this produces a determinate compatible system of linear equa- 
tions: {Z)(^)F(pi) = 0} 

i,h&Hl-> ^"^"^ indeterminate compatible system of linear equations 
{D^h)F{pi) = v,^h} i,h€Hi, where the rank of this system depends on the number of gaps. 

Section four is devoted to exhibit the algorithms we develop and Section five contains some 
examples of how these algorithms work. In particular we are interested in showing that for 
a given real function the addition of conditions on higher order derivatives produces a better 
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approximation. 



1 Multivariate Lagrange interpolation 

The problem of multivariate interpolation consists in given t different elements, pi, . . . ,pt S 
A"(K) in the affine n-dimensional space, and t elements vi, . . . ,vt £ K , determining a poly- 
nomial F G K[Xi, . . . , Xn] such that F{pi) = Vi for every i = 1, . . . , t. 

The following algorithm allows us to determine a polynomial F satisfying this property; 
see [SI [7]. For any index i we consider a linear operator Li : K[Xi, . . . , X„] — > K defined as: 

L,(G) = G{p^), for every G G K[Xi, . . .,Xn]. 

The kernel of Li is Ker(Li) = {G G K[Xi,...,Xn] \ G{pi) = 0}, i.e., Ker(Li) = (Xi - 
Pi^i, . . . , Xn — Pi.n), where pi = {pij, ■ ■ ■ ,Pi,n)- Thus Ker(Lj) is exactly the ideal of the point 
Pi- 

The ideal / = n*^^Ker(Lj) = I{{pi, . . . ,Pt}), is the ideal of the finite set {pi, . . . ,pt}- As a 
consequence the quotient ring K[Xi, . . . ,Xn]/I is a finite dimensional vector space over the 
field K. Therefore there exists a finite basis of K[Xi, . . . , Xn]/L 

There are different methods to compute effectively this basis. Indeed, we consider a monomial 
order in N", hence every non-zero polynomial F G K[Xi, . . . ,Xn] can be written uniquely 
SasAeN" where A is a finite set, ^ ka & K and X°' = X'^^ • • • X^" when- 

ever a = {ai, . . . ,an) G N". We call exp(F) the maximum in A, and define Exp(/) = 
{exp(F) G N" I F G /}. A Groebner basis if I is a set {Gi,...,Gs} C / such that 
Exp(/) = {exp(Gi), . . . , exp(Gs)} + N". The existence of a Groebner basis for each ideal 
I C K[Xi,...,Xn] is well known as it is the uniqueness of a reduced Groebner basis. 
See [4]. The quotient ring K[Xi, . . . ,Xn]/I has finite dimension over K if and only if the 
set B = N"' \ Exp(/) is finite, and it can be taken as an index set for a vector space basis of 
K[Xi, Xn]/L Indeed, the set {X^ + / | /? G 5} is a basis. 

For any polynomial F, which is a solution of the interpolation problem, i.e., F{pi) = Vi for 
any index i = 1, . . . , t, we may write 

F + I = Y^ kpX^ + / for some kp G K. 

Let Fq be the polynomial X^^gg kpX^ , then F{pi) = FQ{pi), for any index i = 1, . . . ,t. Hence 
Fq is an interpolating polynomial which we may compute by solving the linear equation system 
in the unknowns {^/jj/Jeg: 

^ A:/3pf = Vi, i = l,...,t\ 
/3eB I 
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Observe that the uniqueness of the interpolating polynomial is determined modulus the ideal 

1 as if Fi, F2 are interpolating polynomials for any index z = 1, . . . , t we have Fi{pi) = Vi = 
F2{pi), hence F1—F2 S /. We remark that different monomial orders give, in general, different 
interpolating polynomials. 

In order to study the existence of the interpolating polynomial we address to the next sections 
in which we study the multivariate Hermite and Birkhoff interpolation problem. 

2 Multivariate Hermite interpolation 

In addition to the conditions in the Lagrange interpolation problem, it is of interest, some- 
times, to impose some extra conditions in order to put more information in the interpolating 
polynomial. Derivatives is one of the tools we shall use. Intuitively, in the univariate case, 
if we have two different points pi and p2, two values vi and V2, and we want to compute an 
interpolating polynomial F such that F{pi) = Vi we may proceed in the usual way. But if the 
point p2 tends to pi, we would like to use the derivative in pi instead p2- Thus the univariate 
Hermite interpolation deals with the problem of computing an interpolating polynomial F 
such that given different elements pi, . . . ,pt £ K , and elements vq^i, . . . , fo,t, wi^i, . . . , vi^t G K 
the following relations hold: 

n^r=vZ } for every z = l,...,i. 

We may also consider derivatives of higher order. In this case we need to impose some extra 
condition in order to assure we deal with ideals in the polynomial ring. See remarks below. 

Let us describe the multivariate Hermite interpolation problem. 

Let pi,...,pt £ A"(i('), such that pi / pj if i 7^ j, linearly independent unitary vectors 
^ foi' ^ = ji = 1, and elements vi,...,vt G K, Vij.^h G K, 

i = 1, . . . , t, ji = 1, . . . , s(i), h = 1, . . . , s{i,ji). The multivariate Hermite interpolation 
problem consists in determining a polynomial F G K[Xi, . . . ,X„] such that: 

F{pi) = Vi, for every i = 1, . . . ,t; and 

-^[/'l ^(Pi) = '"i,ji,h^ foi' every i = l,...,t, jj = 1, . . . , s{i), h = 1, . . . , s{i,ji) 

In order to determine such a polynomial F let us first consider the problem of determining 
all polynomials F G K[Xi, . . . , X„] such that: 

F{pi) = 0, for every i = 1, . . . ,t; and 

> (2) 

^u,\ ^(Pi) = 0' fo'^ eve'^y i = l,...,t, ji = l,..., s(z), h = l,..., s{i,ji) 
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To study these polynomials we consider a single point, say pi, we assume pi = 0, and the 
system 

F{p^) = 0, 



^uh^(Pi) = 0' fo"^ ^v^'^y ji = 1, • • • , s{i), /i = 1, . . . , s{i,ji) 



(3) 



The sum of two roots of this system also is, and for any root F and any G G K[Xi, . . . , Xn] 
we have: 

{FG){pi) = F{p,)G{pi) = QG{pi) = 0. 

k=0 

Thus the set of all roots of system ([3]) constitutes an ideal, say Ij. We claim the ideal Ij is 
co-finite. In fact, contains the ideal Jj = {X^^ ■ ■ ■ X^" \ ei + • • • + e„ = s{i,ji) + 1). 

Now we may study the set of roots of system ([2]), i.e., the intersection I = IiH . . . H It- Hence 
it is an ideal of K[Xi, . . . ,Xn]. We have that I is co-finite as / is an intersection of finitely 
many co-finite ideals. 

If G is a reduced Groebner basis of / and Exp(/) C is the mono-ideal of all exponents 
of elements in I, a basis of the vector space K[Xi, . . . ,Xn]/I is parameterized by i3 = N" \ 
Exp(/). Indeed, if we consider the basis {X^ \ (3 G B}, a generic element in K[Xi, . . . , X„]// 
has a representative in the form 

^k,X^ 

for some fc^ S K. Every condition in ([T]) produces a linear equation with unknowns in 
{kp I (3 £ B}. Thus we obtain a system of linear equations: 

(E/3eB ^p^^){Pi) = Vu for every z = 1, . . . , t; and 

F>u!^ (E/3GB^/3^^)bi) = '"i,juh, for every i = 1, . . . , t, ji = 1, . . . ,s(i), h = I, . . .,s{i,ji) 

(4) 

Any solution of system gives an interpolating polynomial. 
Observe that if this solution exists, it is unique modulo the ideal /. 



2.1 Generalized multivariate Hermite interpolation 

We may also consider a more general version of the Hermite interpolation problem. Indeed, 
we consider elements pi, . . . ,pt G A"(i^) such that pi / pj if i ^ j, linearly independent 
unitary vectors Uij. G K^, i = l,...,t, ji = l,...,s{i) and sets Hi C p^^W such that if 
h G Hi, and h - e£ W^'^ , then h - e £ Hi for every e G N'*^*). 
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For any h = (hi, . . . , £ Hi and any polynomial F £ K[Xi, . . . , X„] we define 

With this notation, the generalized multivariate Hermite problem consists in, given {vi^h I 
h G Hi} C determining a polynomial F E -^'[^i, • • • , X„] such that: 

D^^^F{pi) = Vi^h, for every i = 1, . . . ,t, h £ Hi } (5) 

Observe that for any index i = l,...,t we may prove, mutatis-mutandi, that the set of 
solutions of 

D^'''^F{pi) = 0, for every i = I, . . . ,t, h £ Hi } (6) 

is an ideal, say Jj. Hence the algorithm may be built as we exemplified above. This will be 
the algorithm we develop in Section HI 

Let us resume in the following theorem these results. 

Theorem 2.1. With the notation in this section there exists a unique interpolating polynomial 
F, modulo the ideal I, such that 

D^^^Fipi) = Vi^h, for every i = l,...,t; h £ Hi^ 

Proof. We only need to prove the existence of such an interpolating polynomial F; it is a 
direct consequence of Corollary (13. 2p . □ 

3 Multivariate BirkhofF interpolation 

The multivariate Hermite interpolation method does not work if there exist some derivative 
gaps, as in this case we can not assure li is an ideal. This problem may be solved if we fill 
these gaps in order to compute the ideal associated to each point pi. Let us perform this 
process and show how to compute this new ideal. 

As in subsection ()2.ip . let us consider points pi, . . . ,pt such that pi / pj if i / j, linearly 
independent unitary vectors Uij- £ K^, i = 1, . . . , t, jj = 1, . . . , s{i) and sets Hi C As 
usual, for any index i, any h £ Hi and any polynomial F £ K[Xi, . . . , Xn] we define 

z)Wf = • • • z?!f^w V. 

The multivariate Birkhoff interpolation problem consists in given {vi^h \ h £ Hi} C K 
determining a polynomial F £ K[Xi, . . . , Xn] such that 

D^^^F{pi) = Vi^h, for every i = I, . . . ,t, h £ Hi} 
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If we fix an index i, we are interested in the set of all polynomials 

{f G K[Xi, ...,Xn]\ D^''^F{p^) = 0, for all h £ 



Contrary to the theory developed in subsection (|2.ip . this set of polynomials is not an ideal. 
To remedy this hitch we proceed as follows. First we define bi the maximum in the set 

s{i) 

\h\ ■.= ^hi\ h = {hi)iGH, 
1=1 

and define Hi = {h & W^^\ such that \h\ < bi}. Second we consider the polynomial equations 

D'^^'^F{pi) = 0, for every /i G ^} (7) 

Let Ji the set of all polynomials satisfying these equations. We obtain that Jj is an ideal of 

K\Xi,...,Xn]. 



Similar arguments, to those developed in section give us that the ideal Ji is cofinite. 
Hence if we define J = Ji H . . . n Jt, then J is a cofinite ideal of K[Xi, . . . , 

Let G be a Groebner basis of J and B the complement in N" of Exp(J), then B has finitely 
many elements and {X^ | /3 G ;B} is a vector space basis of the quotient K[Xi^ . . . ,Xn]/ J- 
So every element in K[Xi, . . . ,X„]/J has a unique representation of the shape X^^gg kpXl^ , 
where G K. 

In order to compute an interpolating polynomial F which is a solution of the multivariate 
Birkhoff interpolating problem we need to solve the linear equation system 

D^^^ ^ k^x" I (k) = Vi^h, i = i,...,t, heHi 

Theorem 3.1. With the notation in this section there exists an interpolating polynomial F 
such that 

D^'^'^F{pi) = Vi^h, for every i = l,...,t; h £ Hi^ 

Proof. For any index i = 1, . . . , t we consider the ideal Jj. This ideal can be computed in an 
easy way. Indeed, if we assume pi = 0, then Jj is the ideal generated by all monomials X^ 
such that |e| > 6j + 1. In particular, a minimal set of generators of Jj is 

{X'\ \e\=bi + l}. 

The codimension of Ji is finite, say di. 

We claim the codimension of J = Ji H ... ri Jt is d := di + ■■■ + df. 
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First we show that Ji + Jj = K[Xi, . . . ,Xn] whenever i / j. Indeed, in this case there is a 
component, say 1, such that Pj,i- Since (^i— G Ji, {Xi—pj^i)^^ G Jj and Xi—pi^i, 

Xi~Pj,i are coprime, then Ji + Jj = K[Xi, . . . ,Xn]- We claim (JiH J2) + J3 = K[Xi, . . . ,Xn]- 
Indeed, there are xi^s £ Ji, ^2,3 G J2 and 2:3^1, ^3^2 G J3 such that xi^3 + X3^i = 1 = ^2,3 + 2:3^2, 
then 1 = a;i,3a;2,3 + (3^1,32^3,2 + a;3,ia;2,3 + 3:3,1x^^2) G (Ji n J2) + -^3- By induction we obtain 
( Ji n . . . n Ji_i) + Ji = K[Xi, ...,Xn] for any 2 < i < t. 

Now the result follows, by induction on i, from the following isomorphisms: 

K[Xi,...,Xn]/{Jin...nJi) ^ K[Xi,...,Xn] 
Ji/{Ji r\...r\Ji) Ji 

Ji _ Ji n . . . n J/_i + J/ _ K[Xi, Xn] 



( Ji n . . . n Ji) Ji n . . . n Ji_i Ji n . . . n Ji„i 

Observe that #{B) = d = di ^ h dt. 

We consider the linear system 

k^X'^^ {pi) = 0, for every i = 1, . . . , i, G (8) 

This system has d equations as it has di equations for every index i. The matrix of this system 
is a square matrix. Since the system has unique solution, this matrix is regular. 

Therefore if we consider the linear system 

^^'^ ^^^'^ ^P'^ = fo'' ^^^^y i = l,...,t, heH)j (9) 

it is an indeterminate compatible linear system. In particular there is at least a solution. 
Hence there is an interpolating polynomial which is a solution to the multivariate Birkhoff 
interpolation problem.. □ 



We remark that the solution to the Birkhoff problem is not unique, in general if there are 
some derivative gaps. 

As a consequence of this Theorem we have the following result which completes the proof of 
Theorem (f^ 

Corollary 3.2. There is an interpolating polynomial solving the multivariate Hermite inter- 
polation problem. 
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3.1 Final comments 



Increasing the number of restrictions 

In the three multivariate interpolation methods we have studied if we add a new point, it is 
not necessary, in order to compute an interpolating polynomial, to start again the algorithm 
from the very beginning. Let us show this situation in the case of the multivariate Birkhoff 
interpolation method. Let pt+i be a new point, different frompi, . . . ,pt, ?7t+i,i, • • • , ?7t+i^s((+i) 
unitary vectors in K", Ht+i C I^'^C'+i) and elements {vt+i^h ^ K \ h € Hf+i}. To compute 
the new interpolating polynomial we only need to compute the following elements: 

(1) the new ideal Jt+i by computing a Groebner basis G^+i; 

(2) the intersection J H Jt+i by computing a Groebner basis G, using the Groebner basis of 
J and Gj_|_i; 

(3) a basis of K[Xi, J n Jt+i) using N" \ Exp(G) and 

(4) solve the linear system 

D^^'^F{p.i) = Vi^h, for every i = 1, . . . , t + 1; h £ Hi^ 

We proceed in the same way when we add some extra restriction on high degree in the 
derivatives at the points pi, . . . ,pt. 

Different orders in N"^ produce a different shape of the interpolating polynomial 

It is of interest to observe that different monomial orders in N" give different shape of the 
interpolation polynomial. For instance, the lexicographic monomial order with the ordering 
of unknowns Xi > X2 > • • • > Xn produces that the interpolating polynomial has greater 
degree in the greatest labeling unknowns. Otherwise, the graded (reverse) lexicographic order 
produces interpolating polynomials in which all unknowns have similar degree. The reason is 
that the shape of the set N" \ Exp(G) strongly depends on the chosen monomial order. 

Complexity of the algorithm 

There are different factors in order to perform a quick algorithm. 

First we need to compute a Groebner basis for each ideal /j or Jj, depending of the methods. 
This process is well established; the speed depends of the kind of numbers we use in the 
coordinates of points. For integer numbers it works in an "acceptable" way and for other kind 
of numbers, rational, reals, complexes, single or multiple precision, etc., it is more and more 
slow. 
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The second step is computing a Groebner basis of / or J. We observe the best method is to 
make it as follows: Ii Ci I2, h ri I2 ri I3, . . . The reason is that anyway we need to perform a 
total of t — 1 intersections of ideals and the described method realizes at the i-th steep the 
intersection (Ii H • • • H /j) H Ij+i, where at least /j+i has a simple description. 

The third step is to determine a general interpolating polynomial using N" \ Exp(G) as set of 
indices, and after that to solve a linear system. It is at the last point where the type of data 
is again determinant in order to get quickly a solution. We also may solve numerically this 
system in order to get an approximation to the final solution. 

Existence and uniqueness of the interpolating polynomial 

As we mentioned earlier different monomial orders in N" produce different interpolation poly- 
nomials, and we also proved in Theorem (13. 1|) the existence of a solution to the multivariate 
Birkhoff interpolation problem. Otherwise the uniqueness of an interpolating polynomial de- 
pends on the restrictions we have, but in any case it can be established only modulo the ideal 
/ or J, depending of the method. In Theorem (12. see also Corollary (|3.2|) . the uniqueness 
was established for the multivariate Hermite interpolation problem. In the Birkhoff interpo- 
lation problem uniqueness fails as we show in the section devoted to examples at the end of 
this paper; the reason is the existence of some derivative gaps at some point. 

4 Multivariate Interpolation Algorithm 

We describe the algorithm for the Hermite interpolation problem. With the same notation 
that in subsection 12.11 let pi,...,pt G A"(K), such that pi 7^ pj if i ^ j, unitary vectors 
Uij. £ K^, for i = 1, ... ,t, ji = 1, ... , s(i), and elements {vi^h I h E Hi, i = 1, . . . ,t} C K . 

The algorithm has three blocks. 

(I) The first one determines a Groebner basis of the ideals Jj of all roots of system ([6]), 

(II) The second one determines a Groebner basis of I = Ii H ... Ci It. 

(III) The third one gives the interpolated polynomial F satisfying conditions in ([5]) . 

We now describe explicitly the developed algorithm for computing a Hermite interpolating 
polynomial. 

4.1 Block I 

Compute a Groebner Basis of every ideal /j, for any index i = 1, . . . ,t. 
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LETp:=pi = (xi,...,x„)GA"(i^) 

LET r := h = max{|/i| := Y.f=i hj \ h = {hj)j £ Hi} 

LET s := s{i) 

LET H := Hi 

LET GEN{I) := 

IF r = THEN 

GEN{I)={X,-X,,...,Xn-Xn} 

ELSE 

FOR / = 1 TO r 

LET a{l) := {a G N" | \a\ := ai -\ \- a„ = 1} 

LETP:=X:H=i^(«)^r---^^ 
LET Deri ■= {h e H \ \h\ = 1} 
LET m := #{Deri) 
FOR A; = 1 TO m 
LET h{k) E Deri 

LET Pk := D^'^^^P = E|a|=/Cfc(a)^(a) 
ENDFOR 

LET / : ^ i^-, /((^(a)),e«(o) = (E|«|=z Cfc(a)^(a))^^, 

COMPUTE a basis of the kernel of the hnear map / 

FOR g = 1 TO dim{Ker{f)) 

LET dq{a) be the g-th element in the basis of Ker{f) 

ADD TO GEN{I) the element dq{a){Xi - • • • (X„ - 

ENDFOR 
ENDFOR 

ADD TO GEN (I) all the monomials of degree r + 1 in Xi — xi, . . . , Xn — Xn 
ENDIF 

COMPUTE a Groebner basis IB of {GEN{I)) 
LET Ij = (B). 



4.2 Block II 

Compute / = /i n . . . n giving a Groebner basis. 
LETI:= {Ii,...,It} 

LET O the graded reverse lexicographic order on X = {Xi, . . . , Xn} 
LET V the lexicographic order on {Xo,X} related to O, 

where Xq is the auxiliary variable to realize the intersection through the elimination method 
LET I = Ii 
FOR k = 2T0t 

LET / =INTERSECTI0N(/,4,J') 
ENDFOR 
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INTERSECTION(/, J, provides the ideal intersection of / and J giving a Groebner basis 
of it with respect to the order V. 

4.3 Block III 

Let us to show the algorithm to find the set N'^\Exp(/). After that we only need to establish 
the system of linear equations in jS]) and solve it to find the interpolating polynomial. 

LET I := {gi, . . . , gk) , where {gi, ■ ■ ■ , gk} is the Groebner basis obtained in Block II 

LET (at,i, . . . , Q!t,n) := The leading exponent of gt 

LET rrii := max{ai^i, . . . , ak^i} 

LET p := (mi - 1, . . . , m„ - 1) 

COMPUTE exp := {/i G N" I p - /i G N"} 

LET b:=0 

FOR i=l TO #(exp) 

LET bt := 

LET j := 

WHILE (i < k AND bt = 0) 
j + + 

IF {expi^s >= ai,s,Vs = l,..,n) THEN 

bt = {expi} 
ENDIF 
ENDWHILE 
b = b\Jbt 
ENDFOR 
LET exp := exp — b 

Therefore we take 

F=Y^ A{a)X». 

aSexp 

Remark 4.1. The Birkhoff interpolation problem admits a similar algorithm in which blocks 
II and III are applied in the same way and for block I we have that, in the notation used in 
this paper, a minimal set of generators of the ideal Jj is directly given by 

{{Xi - pi^iY^ ...{Xn- Pi^nT" I ei + • • • + e„ = 6i + 1}. 



5 Examples 

Example 5.1. Let us show an example in which the use of derivatives of higher order gives 
a more exact approximation to the surface we are studying. We begin with a ellipsoid E of 
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equation 

F{X,Y) = {cos{X)cos{Y),cos{X)sin{Y),'isin{X)) 
and fix eight points in E, for instance 

pi = (2,0), P2 = (2,7r/2), p3 = (2,vr), P4 = (2,37r/2), 

P5 = (3, vr/4), p6 = (3, 37r/4), = (3, 57r/4) and pg = (3, 77r/4). 

We have the following picture: 




We may apply multivariate Lagrange interpolation with the following values 



Pi 


= (2,0), 


Vl 


= (cos(2),0,3sm(2)) 


P2 


= (2,vr/2), 


V2 


= (0,cos(2),3sm(2)) 


P3 


= (2,vr). 


V3 


= (-cos(2),0,2sm(2)) 


P4 


= (2,37r/2). 


■U4 


= (0, -cos(2),3sm(2)) 


P5 


= (3,7r/4), 


V5 


= (cos(3)/\/^,cos(3)/\/2, 3sm(3)) 


P% 


= (3,37r/4), 


V6 


= (-cos(3)/V2,cos(3)/V2, 3sm(3)) 


P7 


= (3,57r/4), 


vr 


= (-cos(3)/ ^2, -cos{3) 3sm(3)) 


P8 


= (3,77r/4), . 


V8 


= (cos(3)/V2, -cos(3)/\/2, 3sm(3)) 



and represent this polynomial. Thus we obtain: 
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Now we try to compute the multivariate Hermite polynomial with derivatives up to order one 
The table of data is: 



Pi 


F{p^) 




D(o,i)F{pi) 


(2,0) 

(2,f) 
(2,^) 
(2,f) 

(34) 

(3,f) 

(3,f) 

(3,^) 


(cos(2),0,3sin(2)) 
(0,cos(2),3sin(2)) 
(-cos(2),0,3sin(2)) 
(0,-cos(2),3sin(2)) 

/ cos(3) cos(3) r, ■ /o\\ 

( ' '3sm(3)j 

/ cos(3) cos{3) o„:„^o\'\ 

( ^°;f),3sin(3)) 

/ cos{3) cos{3) Q„:„/Q\\ 

[ V2 ' V2 '-^sm^djj 


(-sin(2),0,3cos(2)) 

(0,-sin(2),3cos(2)) 

(sin(2),0,3cos(2)) 

(0,sin(2),3cos(2)) 

/ sin|)^ -;;|),3cos(3)) 

(-'^''^'3cos(3)) 


(0,cos(2),0) 
(-cos(2),0,0) 
(0,-cos(2),0) 
(cos(2),0,0) 

/ cos(3) cos(3) n\ 

V V2 ' V2 

1 cos(3) cos(3) (-,\ 

\ ' 72 'V 
/ cos(3) cos(3) n\ 

I V2 ' V2 
I cos(3) cos(3) (-.A 



(10) 



If we represent the Hermite interpolation polynomial we obtain: 
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Observe that in this case the new surface is a better approximation to the original one. 

As a third guest we compute the multivariate Hermite polynomial with derivative up to order 
two. In this case we need to add new columns to the table given in (fTOjl . 



Pi 






Al,0)Ao,l)^(Pi) 


Ao,l)^(0,l)^(Pi) 


(2,0) 

(2,f) 
(2,vr) 

(2,¥) 

(3,f) 

(3,f) 

(3,^) 

(3,^) 




(-cos(2),0,-3sin(2)) 
(0, -cos(2),-3sin(2)) 
(cos(2),0, -3sin(2)) 
(0,cos(2),-3sin(2)) 

/ cos{3) cos{3) q„;„|'oN'\ 

I ' ' -^sm^cijj 

/ cos(3) cos(3) o ■ ro\\ 

( ' ' 3sm(3)j 
/co^3)^cos|)^ 3sin(3)) 

/ cos{3) cos(3) Q„;„,'q^^ 


(0,-sin(2),0) 
(sin(2),0,0) 
(0,sin(2),0) 
(-sin(2),0,0) 

/ sm(3) sm(3) p|\ 
I v/2 ' 72 '^j 
/ sm(3) sin(3) p|\ 
I ' v/2 
/ sin(3) sin(3) 
[ V2 ' V2 '^J 
1 sin{3) sin{3) rC\ 
\ ^ ' V2 


(-cos(2),0,0) 
(0,-cos(2),0) 
(cos(2),0,0) 
(0,cos(2),0) 

/ cos(3) cos{3) p|\ 

\ ' V2 
1 cos(3) cos(3) (-,\ 

/ cos(3) cos(3) n\ 
I V2 ' V2 
1 cos(3) cos(3) (-.A 
I V2 ' V2 



(11) 



If we represent the multivariate Hermite interpolation polynomial we obtain: 
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This representation is a good approximation to the original function. 



Problem 5.2. It would be of interest to measure the error produced by each polynomial inter- 
polation according to the order of the involved derivatives. Nevertheless we have chosen this 
graphical approach to exemplify the goodness of the multivariate Hermite interpolation method 
and postpone a deep study of errors to a forthcoming work. 

Example 5.3. Let us show a new example in which we increase the number of variables, from 
two to four. We also use variables in different types: polynomial, logarithmic and fractional. 
We are interested in studying the following function: 

Observe that this function is of polynomial type in T,Z, of logarithmic type in 1/X and 
of fractional type in X,Y. We propose to study the error produced by the interpolating 
polynomial when we consider a grid of points and evaluate the derivatives at these points up 
to order 2, and show this error through the graphical representation of sections fixing pairs 
of two variables. 

Let us consider the following grid: {{x,y,z,t) G | x,y,z,t £ {1,2,3}} and evaluate the 
functions: 

-f^i -C'{l,0,0,0)-^> o,0)-^> -^^(1, 1,1,0) -C)(l,l,l,l)-^7 

-0(1,0,0,0) -^'(l, 0,0,0) -^5 -C>{1, 1,0,0)-C'(i, 1,0,0)-^; -0(1, 1,0,0)-C'(1, 0,0,0)-^ 

in the grid. We call P the interpolating polynomial. 

In the following table we represent the sections of F and P for the following values of x G 
[l,3],y = 2,zG [l,3],t = 2. 
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F(X,Y,Z,T) 



P(X,Y,Z,T) 



F- P 



In this case the maximum error is smaller than 0.000015, and it is reached next to x = 2.75 
and z = 1 or 3. We have a similar behavior when we fix x and y S [1,3]; in this case 
the maximum error is reached next to y = 1.4. Hence we are interested in fixing these 
two variables in these values. The representation of this section, for the following values: 
X = 2.75, y = 1.4, z, t G [1,3]. 




F(X,Y,Z,T) P(X,Y,Z,T) F-P 



In the case of z,t the maximum error is reached next to z,t = 1 or 3. If we fix these values 
z = 1, t = 3 and allow x,y to vary in [1,3], we obtain the following representation: 
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F(X,Y,Z,T) 



P(X,Y,Z,T) 



F- P 



In this case the maximum error is over 0.001. 

As a final example let us study a little variation of the former one in which we apply the 
Birkhoff interpolating method. 

Example 5.4. Take the real function G{X, Y, Z) = F{X, Y, Z, 2), the grid of points {(x, y, z) \ 
x,y,z G {1,2, 3}} and evaluate the functions: 

G, D(^i i QjG, ^(i^i^i)^, -D(i^o,o)-D(i,o,o)G', -D(i^o,o)-^^(i,i,o)G', D(i^i^o)-^^(i,i,o)G'. 

In order to determine a Birkhoff interpolating polynomial we compute the reduced Groebner 
basis of the ideal J, see Section [Sj this ideal has codimension 27 x (1 + 3 + 6) = 270. Each gap 
in the Birkhoff conditions produces a parameter in the interpolating polynomial; in this case 
we have 27 x 4 = 108 gaps which coincide with the number of parameters. Hence the general 
Birkhoff interpolating polynomial is determined as a solution of a system of 162 independent 
linear equations, whose expression we omit because its large size. 
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